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Abstract — In this paper, a transmission strategy of fountain 
codes over cooperative relay networks is proposed. Wlien more 
tlian one relay nodes are available, we apply network coding 
to fountain-coded packets. By doing this, partial information 
is made available to the destination node about the upcoming 
message block. It is therefore able to reduce the required number 
of transmissions over erasure channels, hence increasing the 
effective throughput. Its application to wireless channels with 
Rayleigh fading and AWGN noise is also analysed, whereby the 
role of analogue network coding and optimal weight selection is 
demonstrated. 

I. Introduction 

Fountain code |[T] and cooperative communication ||2l are 
two transmission strategies which are gaining popularity in 
recent years. In fountain codes, message bits are grouped 
into blocks, each containing several packets. Encoding is 
performed by taking random linear combination of the packets 
within each block over a Galois Field (typically GF{2)). As 
opposed to other fixed rate scheme, with fountain code, the 
source continuously transmits encoded packets until positive 
acknowledgement is received. Hence, its optimality is guaran- 
teed for erasure channels, regardless of the erasure probability. 

Cooperation, on the other hand, improves transmission 
quality by making use of neighbouring nodes to forward the 
message to destination. By creating multiple paths between 
source and destination (each subjected to independent fading), 
a diversity advantage can be exploited. Application of fountain 
codes in cooperative network have also been studied. For 
example, reference (3] and |4| proposed to combine distributed 
space time block coding (DSTBC) and fountain code for 
single carrier and multiple carrier transmission respectively, 
and showed that extra diversity gain can be achieved. In fS^, 
direct application of fountain code in cooperative network 
is analysed. It was shown that careful degree distribution 
design is necessary to ensure good decoding performance. 
Alternatively, using conventional degree distribution, the en- 
coding/decoding process at the relay node can be modified to 
cater for the online re-coding requirement, as discussed in 0. 

Although it has been shown that performance improvement 
can be achieved using fountain codes in cooperative networks 
Q, this advantage brings about extra complexity, especially 
when more than one relay node is involved. Motivated to 
address this issue, recently the authors have proposed an 
amplitude modulation scheme for fountain code transmission 
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Fig. 1 . System diagram of two-relay cooperative erasure network 

over multiple relay cooperative networks fSl. However, the 
scheme is developed for erasure channels, and it is not directly 
applicable to wireless fading channels. The focus of this paper 
is to find an alternative strategy to tackle the above issue. Here, 
a novel scheme that combines network coding L9J and fountain 
codes transmission is proposed. In erasure channel, the scheme 
applies digital network coding onto fountain encoded packets 
of two consecutive blocks, and allows source node to trans- 
mit together with the successful relay. Whereas in wireless 
channel, the scheme employs analogue network coding with 
appropriate power allocation. The performance of the scheme 
is then analysed numerically, and it is shown to improve the 
overall throughput in both types of channel. 

The rest of this paper is organised as follows. Section 
nil describes the system model. The proposed transmission 
scheme and its application into wireless channel are given in 
nni and |IV] respectively. Numerical results are then presented 
in Section |V] Finally, Section |VT] gives concluding remarks. 

II. System Model 

A half duplex cooperative network with one source S, 
one destination D, and two relay nodes (denoted as Ri and 
R2) as depicted in Figure [T] is considerecu. At any one 
time, a node can either transmit or receive, but not both 
simultaneously. The channel between any given two nodes is 
modelled as erasure channel, with a superscript indicating the 
nodes under consideration (e.g., P^^ is used to indicate the 
erasure probability of the channel between S and D). 

Message bits at the source are grouped into blocks of K 
packets, each composed of m information bits. Fountain code 

'Although the discussion presented in this paper is mainly for two relay 
nodes scenario, the results can be extended into general number of relays. 



is then applied to the K packets in which linear combination 
of randomly selected d (generated following some degree 
distribution p{d)) out of K packets are transmitted. For sim- 
plicity, random linear fountain code is considered throughout 
the analysis. Therefore, the degree distribution used to generate 
encoded packets can be expressed as: 
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which can be approximated as p{d) ?» Cf /2^ for large K. 
Here, C™ denotes the number of combinations for selecting 
n out of m elements. The analysis presented in this paper is 
independent of the actual degree distribution used, therefore 
the same technique can be applied to other (more practical) 
degree distributions such as Robust Soliton Distribution 1,10,1 . 

III. Transmission Schemes 

A. Direct Transmission (Without Relay) 

As a baseline comparison, direct transmission of fountain 
code from 5 to Z) is considered. In this case, the number of 
encoded packets received unerased at D (denoted as N) is 
random, and it is related to the number of transmitted packets 
(Af , where < iV < M) through binomial distribution with 
parameter P^^ as follows: 

Bm.PSo (N) = C^^ (1 - Pf ^)^ (^pSD^M-N (2) 

Given that N encoded packets are available at D, the prob- 
ability that the corresponding K x N generator matrix is 
full rank can be calculated. Following the assumption that 
random linear fountain code is used, all 2^^ possible binary 
generator matrices are equiprobableB Hence, the probability 
of successful decoding can be calculated as: 
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The above equation can be interpreted as the cumulative 
distribution function (CDF) of successful decoding after N 
packets are received. The respective probability distribution 
function (PDF) can be expressed as: 

f{N)=F{N)-FiN-l) 







N < K 
N = K (4) 
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Noting that the transmission stops as soon as D is able to 
decode, the PDF of successful decoding after A/ transmissions 
is given as: 



pW(A/,Pf^) 
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where fi is the normalising constant to satisfy the unit-sum 
constraint T, m P^^H^ , P^ " ) = 1, while B{-) and /(•) are 
given in dU and (|4|i respectively. 

^Note that we have included the probability of generating degree column 
to simplify the expression. This approximation is tight for large K. 
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Fig. 2. Illustration of two-phase transmission using naive relaying scheme 
(assuming iJi is the successful relay). 



B. Naive Relaying Scheme 

In the presence of relay node, a natural way to perform 
transmission is through multiple hops using two-phase scheme 
as illustrated in Figure |2l During the first phase, S broadcasts 
fountain encoded packets to all i?'s, as well as D. Since P^^ 
(assumed to be equal for all relays) is smaller than Pf^ due 
to geographical proximity, it is highly likely that one of the 
i?'s is able to decode the message before D. Denoting the 
total number of relays as TZ, the probability that any of them 
is able to decode after j transmissions can be calculated as: 



qU) = Y.cf [p^'\j.p!")\ [Y.p^'H^.p!"T- 
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The above equation holds for any value of TZ, and it determines 
the duration of phase one transmission. Once any of the i?'s 
decodes, the transmission enters the second phase, whereby the 
successful R transmits fountain encoded packets to D. Since 
the last packet from R to D must not be erased, the PDF of 
successful decoding after M transmissions is given as: 



p(2)(Af)=p(i)(M, 
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where the first term represents the case when D is able 
to decode before any of the P's, hence the probability of 
successful decoding falls back to direct transmission case; and 
the second term represents the case when one of the R\ is able 
to decode before Z?, hence the received packets come from 
both S and R during the first and second phase respectively. 
Remark: Due to broadcast nature of the transmission, the 
other relays can receive the packets from the transmitting relay 
during phase two. However, those packets would not be useful, 
as they only contain information about the current message 
block. The moment D is able to decode, it will send positive 
acknowledgement (ACK), and a new protocol cycle starts. 
Here, ideal ACK channel is assumed available for simplicity. 

C. Proposed Network Coded Scheme 

The idea of this scheme is to allow S transmitting partial 
information about the next message block during the sec- 
ond phase, hence enabling the receiving relays to help the 
transmission of subsequent message blocks. The first phase 
transmission in this case is identical to the naive relaying 
scheme. On the second phase, in addition to the successful 
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Fig. 3. Illustration of two-phase transmission using the proposed network 
coded scheme (assuming iJi is the successful relay). 

TABLE I 
Buffer assignment table at destination node during phase two. 
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R, S also transmits the network coded version of fountain 
codes for the current and the next message block as illustrated 
in Figure [5] Here, the fountain code for the current message 
block at the successful R and S are identical, while different 
code can be used for the next message block at S. 

Without loss of generality, let Ri be the successful relay, 
and there are only two relays in the network. Since the channel 
is an ideal erasure channel, the effective packet received at 
i?2 and D is the combination (bit-wise exclusive-or) of the 
packets transmitted from S and i?i. Depending on which of 
those packets are erased, D would store the resulting packet 
to one of the three buffers it maintains as illustrated in Table U 
When the packet from R is the only one unerased, the received 
packet contains information about the current block, and it is 
stored in buffer 1. When the packet from S is the only one 
unerased, the packet contains a network coded version of the 
current and next block, and it is stored in buffer 2. Finally, 
when both packets are unerased, the received packet contains 
information about the next block, and it is stored in buffer 3. 

To enable the receiving nodes to identify which one of the 
three cases has happened, a modification is made to the packet 
structure, as depicted in Figure |4] 



h, hi 



Encoded packet 



flbits mbits 

Fig. 4. Modified packet structure with multiple headers. 

By assigning header hi exclusively to S and /i2 to R, the first 
H bits of the packet can be used to identify which packets has 
been erased. Meanwhile, the remaining m bits are shared by 
both S and R. As long as m >> fi, the overhead caused by 
the additional header information is negligible. 

Comparing this scheme with the naive relaying, several 
observations can be made. Firstly, the useful fraction of 
transmissions during the second phase can be calculated by 
adding the probabilities of the three cases, which is equal to 
(1-P^^^)+(1-Pf ^)P,^^. This is lai-ger than that of the naive 
relaying method, which is only (l—P^^). Secondly, the pack- 
ets received at R2 during phase two contain information about 



the next message block, which can accelerate its transmission 
to D once the current message is decoded. 

To perform the combinatorial analysis of this scheme, the 
following auxiliary function is defined: 






Ef=o'f3M-l.pA^)fnAi + '^) M>0 



where 5I2 is the normalisation factor required to make 

J2ai 9np.pAM) = 1. Function /„p(.) is defined as: 



fnJN) 



F{np) N ^Q 

F{N + np)-F{N + np-l) N>0 



(9) 



which can be evaluated by substitution with equation (|3]l. 
The above auxiliary function g„ p^ (A/) can be interpreted as 
the probability of successful decoding (in point-to-point case) 
after M transmissions, given that D already has Up packets 
available. Denote the number of packets available at i?i, R2, 
and D as rii, n2, and 71,3 respectively. The probability that 
either i?i or R2 is able to decode after j transmissions is: 

(10) 
Correspondingly, the probability that D is able to decode after 
M transmissions can be calculated as: 
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where PJ^* = 1 — {P^^{1 — P^^)) is the equivalent erasure 
probability of the current message block. 

Initially, all ni, 712, and n^ are equal to zero. Whenever any 
of the P's is able to decode before D (let Pi be the successful 
relay), there will be some network coded packets transmitted, 
therefore some information about the next message block will 
be available at D as well as P2. Let j and AI {M > j) 
be the time at which Pi and D decodes the current block 
respectively. Since D can decode the current message block, it 
can remove its contribution from the network coded packets; 
therefore the content of buffer 2 and 3 can be combined to 
provide the information about the next message block. 

During the M — j transmissions in phase two, D would 
have filled buffer 2 or 3 in 1 — P^^ of the time. Number 
of packets 713 for the next message block is then binomially 
distributed according to I3f,[_j psn{nz). Meanwhile, there are 
two scenarios to be considered for P2- When it is also able 
to decode the current block, which occurs with probability: 

i M-j 

7 = X! X! /5j,n^'^(*)/^M-i,Pe='' (^)/»2 {S + t) 

s=0 t=0 

where P^^^ = 1~{P^'^{1- P^^)) is the equivalent erasure 
probability, the content of buffer 2 and 3 can be combined. 



Therefore, the equivalent erasure probability for the next 
message block at R2 during M — j phase two transmissions 
is P^^ ■ Correspondingly, ri2 for the next message block is 
binomially distributed according to l3j^,[_jpsR{n2)- 

For the case when R2 is not able to decode the current 
message block (which happens with probability 1— 7), only the 
packets stored in buffer 3 can be used to provide information 
about the next message block. In this scenario, since only 
(1 - Pe^){l - P^^) = 1 - -Pg''^ of the transmissions go to 
buffer 3, the number of packets 112 for the next message block 
is binomially distributed according to I3j^j_, p<iq3{n2)- 

As far as i?i is concerned, due to half duplex constraint, 
it will not be receiving any packets about the next message 
block, therefore rii = for the next block transmission. Using 
rii, 712, and n^, transmission will then continue with phase one 
transmission again until all message blocks are sent. 

IV. Application to Wireless Channels 

In wireless channel, instead of erasure, fading and AWGN 
(Additive White Gaussian Noise) contribute to channel impair- 
ments. The channel model between S and D in this case can 
be expressed using its discrete-time equivalent model: 



Vd = hsDXs + V 



(12) 



where hso is the Rayleigh flat fading channel coefficient, 
distributed according to CA/'(0, 1/Xsd) (^sd is the path loss 
parameter), ry is the AWGN noise, normally distributed with 
variance a^. xs and yn are the transmitted and received signal 
respectively. For other node pairs, the subscript notation can 
be changed accordingly. Two approaches to extend the earlier 
results into wireless channel scenario are given as follows. 

A. Approach 1: Convert wireless channel into erasure channel 

Here, the idea is to use a Gaussian codeboolu, and map 
every possible 2™+'' encoded packet into one of the valid 
codewords. The encoding and decoding process can be ex- 
plained using standard random coding argument as follows: 

Encoding: Generate 2"^ length-n codewords, comprising 
of typical sequences drawn according to A/'(0, cr^^,). Here, a^^ 
is the available power of the transmitting nodes, which are 
assumed to be equal for all nodes. Then, a one-to-one mapping 
is performed to associate each one of the 2™+'' encoded 
packets to a valid codeword. Transmission of a particular 
encoded packet is done by transmitting the corresponding 
codeword to the channel. 

Decoding: Perform coherent detection (equivalent to divi- 
sion by channel fading h). Decoding is then performed by 
searching in the Gaussian codebook, the length-n sequence 
which is jointly typical with the equalised received sequence. 
The desired encoded packet can then be obtained by inverse 
mapping, using the same table as that in the encoder. 

Following Shannon's theorem of reliable communication, 
since x is drawn from Gaussian codebook and ?/ is normally 

^^The use of Gaussian codebook is to simplify the exposition on how 
the proposed scheme works in wireless channel. In practical scenario, good 
channel code with reasonable encoding/decoding complexity needs to be used. 



distributed (independent to x), the capacity of the above 
channel is C = 0.51og2(l + \hsD\'^(^'ix/'^n)- ^s long as 
R < C, any codeword can be decoded with vanishing error 
probability as the codeword length n grows. This forms the 
first condition for reliable packet transmission. 

The second condition is related to the mapping table used 
to associate each encoded packet with a valid codeword. Since 
the mapping is required to be one-to-one, it is necessary that 
nR > 771+ ij,. Combining with the first condition, the constraint 
for successful packet transmission can be written as: 



0.51og2(l + \hsD\''<JtJ<jl) > {7n + fi)/7i 



(13) 



Whenever this constraint is violated, there is no rate R 
which satisfies both of the above conditions, hence it is not 
guaranteed that the packet can be decoded without error. 

Considering the non-decodable packet as erased, the channel 
can be treated as erasure channel with erasure probability: 



Pf^=Pr [0.51og2(l + \hsD\^<jl/al) < 
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With the above strategy, each point-to-point channel has been 
converted to its equivalent erasure channel. Therefore, the 
earlier analysis on phase one transmission can be readily 
applied. For the phase two transmission, however, both S 
and successful R are transmitting simultaneously; therefore a 
form of multiuser detection is required at D. Since Gaussian 
codebook is used, decoding can be performed successively 
yjj- First, the codeword of one source is decoded. Then it is 
subtracted out from the received signal, followed by decoding 
of the other source. The erasure probabilities in this case is 
determined by the decoding order used. When the channel 
from i? to D is stronger, i.e. |ft.i?Dp > |/iS£>P (which happens 
with probability Xsd/{Xsd + Xjio)), the packet from R is 
decoded first. The corresponding erasure probability is: 



pRD(a) ^ Pj. 
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Xsd/x 
Xrd + Xsd/x 



\hRD?CF^tx 



-) < (m + /i)/n 



\hsD?CTtx + ^'; 



exp (-AflD X CT^/c^te) (15) 



pSDia) ^ p^ [0.51og2(l + \hsD?(jl/al) < (m + M)/n] 

= 1 - exp (-AsD X (yl/<ylx) (16) 

where x = 22('"+^)/" - 1. On the other hand, when |/iflD p < 
\hsD\^, packet from S is decoded first, and correspondingly: 

Xrd/x 



pSD{b) 



XsD + Xjid/x 
Pf^W = l-exp(-A«z.X->L) 



exp (-AsD X fT^/^L) (17) 



(18) 



It is apparent that there is a dependency to the above erasure 
probability. Namely, when packet from R is decoded first, 
erasure probability for packet from S is equal to Pe 
only when packet from R can be successfully decoded (with 
probabiUty 1 — Pe ), otherwise it will not be decodable. 

It should also be noted that when both packets from S and 



TABLE II 

Buffer assignment table at destination during phase two. 
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Fig. 5. Capacity region of network coded cooperative scheme 



R can be decoded, apart from the encoded message of the 
next block, the encoded message of the current block will 
also be decodable. The buffer assignment table as well as 
the equivalent probabilities are given in Table HI] The overall 
performance can then be analysed using the same approach as 
that in ideal erasure channel. 

B. Approach 2: Assume that fountain code achieves capacity 

In this approach, we focus more on the amount of infor- 
mation to be sent across to destination, which is equal to 
Km bits in each block. Denoting the capacity between S and 
D as CsD = 0.51og2 (l + T\hsD\'^^'ixl^n)^ the duration of 
phase one transmission can be calculated as Tgi = Km/CsR- 
The constant T is the SNR-gap used to cater for the im- 
perfection of the actual modulation used as compared to 
Gaussian codewords IIT2l - |fT3l . During this Tsi period, D also 
received CsdTsi bits of information. Therefore, the remaining 
information about the current block that D requires during 
phase two is Km — [CsdTsi) bits. 

Following the proposed scheme, both S and i?i transmit to 
D in phase two. Here, there are several possibilities. In one 
extreme, S can transmit information about the same message 
block (denoted as Xi). The received signal at D is then: 



y = hsDXi + hfiDXi + rj 



(19) 



After coherent detection is performed with respect to the 
combined channel, the resulting capacity can be expressed as: 

C = 0.51og2 (1 + T\hsD + hRD?a''tJcTl) (20) 



With this approach, the duration of phase two transmission is 
minimised. However, no information about the next message 
block is sent, hence similar problem arises. In the other 
extreme, S can transmit information about the next message 
block. The received signal model at the destination is then: 



y = hsDXi+i + hRDXi + r] 



(21) 



It is observed that the above model is a multiple access channel 
model. The capacity region of this channel is shown by the red 
line in Figure|5] As compared to the previous case (S transmits 
only block i information, whose capacity region is reflected as 
the blue curve in Figure |5]l, the rate of transmission about the 
current block is smaller Therefore the duration of phase two 
is also longer. On the other hand, some information about the 
next message block will also get transmitted to both R and D, 
which will expedite the transmission of subsequent message 
blocks. However, the advantage of having some information 
about the next message block versus the penalty due to longer 
phase two duration needs to be studied. 

Having considered the two extreme cases, it is natural to 
analyse the condition between those two extremes. In this case, 
5" transmits a combination of the message from current block 
i and the next block z -t- 1, resulting in the following: 

y = hsD{aXi+ I3xi+i) + hRDXi + ri 

= {ahsD + hRD)xi+ PhsDXi+i+T] (22) 

This scheme is also known as analogue network coding 
or superposition coding ||T41 . To satisfy the transmit power 
limitation, additional constraint a^ + /?^ = 1 is imposed. It 
is apparent that the above channel is another multiple access 
channel, with channel gains {ahso + ^/jd) and {Phso), 
which are used to transmit information about current and the 
next block respectively. The capacity region of this scheme is 
shown by the black curve in Figure |5j 

Several observations can be made about the new capacity 
region. Firstly, as the value of a is increased, the rate for 
transmitting the next message block z + 1 decreases. Secondly, 
as a increases, the rate of transmitting the current message 
block i is increased, reducing the time required in phase two 
transmission. Lastly, another interesting observation is that the 
sum rate capacity (indicated by the doted green line in Figure 
|5]l increases with parameter a, which suggests that larger a 
is better. However, a should not be too large, as otherwise no 
information about the next block can be delivered. 

Taking the union of all capacity regions as parameter a is 
varied, the entire capacity region of network coded transmis- 
sion scheme can be found (depicted as the pink curve in the 
Figure|5]l). Optimal operating point lies in the boundary line of 
this region. However, further investigation is required to find 
the best trade-off point. One good operating point is to set the 
transmission rate for the current message block to be equal 
to the rate when only R transmits during phase two. In this 
way, the duration of phase two transmission is left unchanged. 
Meanwhile, by operating at the boundary line of the capacity 
region, the amount of information about the next message 
block at the destination (and receiving relay) is maximised. 
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Fig. 7. Average number of transmissions for successful decoding over 
wireless channel as a function of SNR 



V. Numerical Analysis 

In this section, the performance of various schemes in terms 
of required number of transmissions are analysed. Due to 
the nature of fountain code, no analysis on bit error rate is 
given, as transmission is allowed to continue until successful 
decoding is achieved at destination. A system with single 
destination and two relay nodes is considered. The erasure 
probability of the direct link is set to P^^=QA, while that of 
the indirect Hnks is set to P^"=P^^=P^^=Q.2. The number 
of packets in each message block is K=\QQ, and transmissions 
are normalised with respect to packet size m. Figure |6] shows 
the histogram as well as the corresponding PDF (calculated 
using equation (fTTt ) of the required number of transmissions 
using the proposed network coded scheme. The PDF of the 
number of transmissions under direct transmission and naive 
relaying schemes (calculated according to equation (|5]l and (jT) 
respectively) are also shown for comparison. 

It is observed from the figure that the proposed scheme 
outperforms both direct transmission and naive relaying, which 
is reflected by the smaller mean in the required transmissions. 
It is to be noted, however, that the variance is larger; which 
signifies that it is more suitable for non delay-constrained 
applications. Another interesting observation is that unlike 
the other two schemes, there is a positive probability that 
successful decoding is possible in less than K transmissions. 
This is attributed from the extra information in the form of 
encoded packets, which are received during the previous block. 

In Figure [T] the average required number of transmissions 
over wireless channel is depicted for the three schemes under 
consideration. Here, the path loss parameter is related to 
distance according to A oc d", where path loss exponent a 
is set to 3, and the source-destination distance {dso) is set to 
20 meters. Meanwhile, the distance of the relay with respect to 
both source and destination is set to dsR = duo = 10.3 meters, 
and the inter-relay distance is dj^jj = 5 meters. From Figure 
|7] it is observed that the schemes exhibit similar behaviour in 
wireless channel (for both approaches), whereby the proposed 
network coded scheme outperforms the other two schemes. 
The difference is more prominent in low SNR regime. In high 
SNR regime, the effective erasure probability is very small, 
hence the performance of the three schemes are comparable. 



VI. Conclusion 

In this paper, a novel transmission scheme of fountain code 
over cooperative relay networks with multiple relay nodes is 
proposed. The scheme combines network coding concept to 
allow source node transmitting partial information about the 
next message block during the current block transmission. This 
extra information will expedite the transmission of subsequent 
blocks, thereby improving the overall throughput. Two possi- 
ble approaches to implement the scheme into wireless channel 
are discussed. It is then shown via numerical analysis that in 
both approaches, the proposed scheme is able to outperform 
the naive relaying scheme, especially in low SNR regime. 
Future work includes a search for optimal operating point as 
well as a development of robust scheme, which could adapt 
the relaying strategy according to channel conditions. 
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